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Abstract. We study the 3x3 elliptic systems V(a(.v)V X u) - V(f)(.Y)V ■ ii) = f, where 
the coefficients a(x) and h(x) are positive scalar functions that are measurable and bounded 
away from zero and infinity. We prove that weak solutions of the above system are Holder 
continuous under some minimal conditions on the inhomogeneous term /. We also present 
some applications and discuss several related topics including estimates of the Green's 
functions and the heat kernels of the above systems. 



1. Introduction 
In this article, we are concerned with the system of equations 

(1.1) V X (a(x)V X m) - V(Z7(x)V ■ m) = / in Q, 

where the unknown u - {u^,u^, w') and the inhomogeneous term / - {f\f^,f^) are vector 
valued functions defined on a (possibly unbounded) domain f2 c R^*, and the coefficients 
a{x) and b{x) are positive scalar functions on Q that are measurable and bounded away 
from zero and infinity. It should be noted from the beginning that the above system dl.lb 
is elliptic. As a matter of fact, the following vector identity 

(1.2) Vx(Vxm)- V(V-m) = -Ah 

implies that in the case when a and b are constants, the above system reduces to 

-aAu + (fl - fo)V(V ■ h) = / in O, 

which (under the assumption that a > and b > 4a/3) becomes the Lame system of 
linearized elastostatics in dimension three; see e.g., Dahlberg et al. Jl]. A special case of 
the system dl.lb is the following system 

(1.3) V X (a(A)V X m) = 0, V-M=0 in Q, 

which arises from Maxwell's equations in a quasi-static electromagnetic field, where the 
displacement of the electric current is neglected; see e.g.. Landau et al. 1211 Ch. VII]. 
In ST4\ . the authors proved that weak solutions of the system ( 11.3b are Holder continuous 
in Q; see also Yin 125). It is an interesting result because in general, weak solutions of 
elliptic systems with bounded measurable coefficients in dimension three or higher are not 
necessarily continuous; see De Giorgi ||6l. Another motivation for studying the system 
dl.lb comes from an interesting article by Giaquinta and Hong IfTOl . where they considered 
the following equations involving differential forms: 

(1.4) d*{o-(x)dA)^0, -£/*A = in Q, 
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where cr{x) € L°°(D.) is a function with cri < cr(x) < 0-2, cri and 0-2 being two positive 
constants, A is a one-form, dA is its exterior differential, and d* denotes the adjoint of 
d (i.e., d* - S, the codifferential). Related to the well-known result of De Giorgi ||5l on 
elliptic equations, they raised an interesting question of whether any weak solution A of the 
equations ( 11.41 ) is Holder continuous in Q. In the three dimensional setting, the equations 
(II.4I 1 becomes the system (11.3b . and thus, in dimension three, a positive answer was given 
in Il4l . Conversely, in terms of differential forms, the system ( 11.11 ) with f - becomes 

d*{a(x)dA) + d(b{x)d*A) = in Q; A = u^dx^ + u^dx^ + w'dx^. 

Similar to the question raised by Giaquinta and Hong IfTOl . it is natural to ask whether weak 
solutions of the above equations are Holder continuous in Q. We hereby thank Marius 
Mitrea for suggesting this question to us. 

In this article, we prove that weak solutions of the system ( II. lb are Holder continuous 
in Q assuming a minimal condition on /, and thus give a positive answer to the above 
question in dimension three; see Theorem 13 . 1 1 below for the precise statement. With this 
Holder estimate at hand, we are able to show that there exists a unique Green's function 
G{x,y) of the system (II. lb in an arbitrary domain O c R^, and it has the natural bound 

\G{x,y)\<N\x-y\-' 

for all x,y e ^ such that Q <\x - y\ < dx A c/,., where dx '■- dist(jic, dQ), a A b := min(a, b), 
and is a constant independent of Q. In particular, when Q - E?, the above estimate 
holds for all x y; see Theorem l6.2l and Remark \6J\ below. It also follows that the heat 
kernel Ki{x,y) of the system ( II. lb exists in any domain Q., and in the case when Q = R^, 
we have the following usual Gaussian bound for Kf(x,yy, see Theorem I? . 1 3 I below : 

\K,(x,y)\ < Nr^'- exp{-K\x - yf/t}, Vr > 0, x,yeR\ 

Another goal of this article is to estabhsh a global Holder estimate for weak solutions of 
the system ( 11.3b in bounded Lipschitz domains. More precisely, we consider the following 
Dirichlet problem 

( V X (fl(x)V xu)=f + Vxg in Q, 
(1.5) i V m = /! in Q, 

y M = on (9Q, 

where Q is a bounded, simply connected Lipschitz domain. We prove that the weak solu- 
tion H of the above problem (11.5b is uniformly Holder continuous in Q. under some suitable 
conditions on the inhomogeneous terms /, g, and h; see Theorem l3.4l for the details. This 
question of global Holder regularity for weak solutions of the system ( 11.31 ) turned out to 
be a rather delicate problem and was not discussed at all in lfT4l . Yin addressed this issue 
in ||25l , but it appears that there is a serious flaw in his proof; he also considered a sim- 
ilar problem with a more general boundary condition in 1261 , but it seems to us that his 
argument there regarding estimate near the boundary has a gap too. Utilizing the above 
mentioned global Holder estimate for weak solutions of the system ( 11.5b . we show that the 
Green's function G(x, y) of the system (11.31 ) in Q has the following global bound: 

\G{x,y)\ < N[dx A\x-y\f[dy Alx-ylflx-yr^---", Vx,yeD., x + y, 

where < a < 1; see Theorem |6.6| for the details. In that case, we also have the following 
global estimate for the heat kernel K,{x, y) of the system ( |1.3b in Q: For all T > 0, there 
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exists a constant such that for all x,y e D. and < f < T, we have 

lA— ■ lA^^^ ■ r^'/'expl-^lx-yp/f), 

VfV|x-y|/ \ ^ftW\x-y\j 

where k > and a e (0, 1) are constants independent of T, and we used the notation 
ay b - max(a, b); see Theorem I?. ISI below. At the moment, it is not clear to us whether or 
not any global Holder estimate is available for weak solutions of the full system (II. lb with 
zero Dirichlet boundary data. 

The organization of the paper is as follows. In Section |2l we introduce some related 
notation and definitions. In Section[3] we state our main theorems and give a few remarks 
concerning extensions of them. The proofs of our main results are given in Section |4] and 
some applications of them are presented in Section |5] We devote Section |6] entirely to the 
study of the Green's functions of the system (II. lb . and Section|7]to the investigation of the 
parabolic system and the heat kernels associated to the system dl.lb . 

2. Notation and Definitions 

2.1. Basic notation. The basic notation used in this article are those employed in Gilbarg 
and Trudinger IfTTII . A Function in bold symbol such as u means that it is a three dimen- 
sional vector-valued function; V ■ u denotes div «, V x « denotes curl u, and Vh denotes the 
gradient matrix of u. Throughout the article, Q denotes a (possibly unbounded) domain in 
'B? (i.e., an open connected set in R^) and 5Q denotes its boundary. For a domain Q. with 
C' boundary dQ., we denote by n the unit outward normal to cJQ. Let L be the operator of 
the form 

Lu := V X {a{x)V x «) - V{b{x)V ■ u) 
whose coefficients are measurable functions on Q satisfying the following condition: 
(2.1) V <a{x), b{x) <v'\ VxeQ, for some v 6 (0, 1 ] . 

For X eQ. and r > 0, we denote B, (x) the open ball of radius r centered at x and 

Q,(x) := Q n Br{x)\ {_dQ.),ix) dO. n Br{x). 
We write S ' cc S if S' has a compact closure in 5 ; 5 ' is strictly contained in S . 

2.2. Function spaces. The Holder spaces C*'"(Q) (C* "(Q)) are defined as the subspaces 
of C*(Q) (C*(Q)) consisting of functions whose ^-th order partial derivatives are uniformly 
Holder continuous (locally Holder continuous) with exponent a in Q. For simplicity we 
write 

C°-"(fi) = C"(n), C"'"(n) = C"(Q), 
with the understanding < ff < 1 whenever this notation is used. We set 

II II II r I ,11 \u{x)-u(y)\ 
m\cc(a) = MQ,a;a = [u\a;Q. + Mq-q. := sup — — + sup|m|. 

x,yen \x y\ n 

For p> 1, we let Lf{Q.) denote the classical Banach space consisting of measurable func- 
tions on Q. that are p-integrable. The norm in L''(Q) is defined by 

\\u\\p;n ^ \\u\\D'(Ci) ^ {ul^dx 

For p > 1 and k a non-negative integer, we let W'^ '^(Q) the usual Sobolev space; i.e. 
W*'^(Q) = {m e LP(n) : D"u e L''{Q) for all lal < k). 
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We denote by Cq(Q.) the set of all functions in C°°(il) with compact support in Q. Some 
other notations are borrowed from Galdi ||8l and Maly and Ziemer ||231 . Setting 

D = = {« e C,7(Q) : V ■ H = in n), 

for q e [1, oo) we denote by Hq{Q) the completion of £)(Q) in the norm of W. The space 
y''^(Q) is defined as the family of all weakly differentiable functions u 6 L^{Q.), whose 
weak derivatives are functions in L^(Q). The space F''^(Q) is endowed with the norm 

:= l|M|lL<>(n) + l|VM||£2(n,. 

If \Q.\ < oo, then Holder's inequality impHes that ^' ^(Q) c W''2(Q). We define Y^^'^iQ.) as 
the closure of C^(Q) in Y^-^{D.). In the case D. = R^ we have Y^-^(R^) = Fo'^CR^). Notice 
that by the Sobolev inequality, it follows that 

(2.2) \\u\y^n) < A^IIVMllL^(n), Vm e Y'^'\n). 

Therefore, we have Wo'^(Q) c Y^/iD.) and Wo'^(Q) = Y^^'^iD.) if |Q| < oo; see ||23l §1.3.4]. 
In particular, if Q is a bounded domain, then we have Y^^'^iQ.) - Wq^{Q.). 

2.3. Lipschitz domain. We say that Q c is a (bounded) Lipschitz domain if 

i) Q is a bounded domain; i.e. 

diamQ :- sup{|x - y\ : x,y e Q) < oo, 

ii) There are constants M and ro > 0, called Lipschitz character of 5Q, such that for each 
P e do., there exists a rigid transformation of coordinates such that P = and 

Q n = = (x', X3) €R^:x3> <p(x')} n B,„; B,„ = B,„(0), 

where : R^ — » R is a Lipschitz function such that ^(0) - 0, with Lipschitz constant 
less than or equal to M; i.e. 

\(p(x') - (p(y')\ < M\x' - y'\, Wx',y' 6 R^. 

2.4. Weak solutions. We say that m is a weak solution in Y^'^(Q.) of the system dl.ll l if 

(2.3) r fl(V X h) ■ (V X <^) + b{W ■ h)(V ■4>)= f f(l>, V<^ e C~(Q). 
Jn Jn 

1 

We say that a function u is a weak solution in Yq''{Q.) of the problem 

V X {a{x)V X h) - V{b(x)V ■ u) ^ f in Q, 
H = on do., 

if u belongs to ^^'^(Q) and satisfies the identity ( 12.31 ). By a weak solution in Fg'^(Q) of the 
problem ( 11.51 ). we mean a function u 6 FjJ'^(Q) satisfying 

(2.5) r fl(VxH)(Vx^)= r f,p + g{Vx4,), V^€Co"(Q) 
Jn Jn 

(2.6) 1 M ■ Vi/r = - 1 ViA € CoiO). 
Jn Jn 

By using the standard elliptic theory, one can easily prove the existence and uniqueness 
of a weak solution of the problem ( |2.4t in Y^^"'(D.) provided / € L^^^iD.). Similarly, if 
/ € //6/5(Q) and g 6 L^(Q.), one can show that there exists a weak solution in ^^■'^(Q) of 
the problem ( IL5b when h - 0; in the more general case when h e L'^'HQ.) and J^h = 0, one 
can show that there exists a unique weak solution in 7^ -(Q) of the problem jl.Si provided 
that Q is a bounded Lipschitz domain; see Appendix for the proofs. 



(2.4) 



ELLIPTIC SYSTEMS WITH MEASURABLE COEFFICIENTS 



5 



3. Main Results 

Our first theorem says that if / e V{Q) with q > 3/2, then weak solutions of the system 
(II. lb are locally Holder continuous in Q. 

Theorem 3.1. Let Q. be a (possibly unbounded) domain in M?. Assume that a{x) and b{x) 
are measurable functions on Q. satisfying the condition ( 12.1b , and that u e Y^'^{D) is a 
weak solution of the system ( II. lb , where f e L'^(Q.) with q > 3/2. Then u is Holder 
continuous in Q, and for all Br — Br(xq) CC Q., we have the following estimate for u: 

(3.2) R"[u]a;B,^, + |H|0;B„,3 < N {R-"^\\u\\i2^^b,) + ^'"'^"ll/llLnB^)) ^ 

where a - a{v, q) € (0, 1) and N - N(v, q) > 0. 

In order to establish a global Holder estimate for weak solutions of the problem ( |1.5b . 
we need to impose some conditions on Q.. We shall assume that is a bounded Lipschitz 
domain whose first homology group Hi{D.; R) is trivial; i.e., 

(3.3) //i(n;R) = 0. 

For example, if Q. is simply connected, then it satisfies the above condition. As mentioned 
in ^2.41 the existence and uniqueness of a weak solution in W^'^(Q) of the problem (11.5b is 
established by a standard argument; see Appendix. 

Theorem 3.4. Let D. c be a bounded Lipschitz domain satisfying the condition ( 13. 3b . 
Let a(x) be a measurable function on Q. satisfying the condition ( 12.1b and u e Wq^{Q.) be 
the weak solution of the problem (11.5b , where f e //^/2(Q), g,h e U'(D)for some q > 3, 
and J^h — 0. Then, u is uniformly Holder continuous in O and satisfies the following 
estimate: 

(3.5) WuWc-m ^ Ni^\f\\u,n-(a) + WgWmsi) + \\h\\mn>) , 

where a - a{v, q, Q) e (0, 1) and N - N{v, q, Q.) > 0. 

Related to the above theorems, several remarks are in order 

Remark 3.6. In Theorem l3.1l one may assume that a{x) is not a scalar function but a 3 x 3 
symmetric matrix valued function satisfying 

v\^\^ <fa(x)^ <v-^\^\^, V^eR^ Vx e Q, for some v e (0, 1]. 

There is no essential change in the proof; see IfTTl . As a matter of fact, one may drop the 
symmetry assumption on a{x) if one assume further that a e L°°(f2). 

Remark 3.7. In Theorem 13. II instead of assuming / e L''{Q.), one may assume that / 
belongs to the Morrey space L^ '* with p = 6/5 and A = 6(1 + 26)/5 for some 6 e (0, 1); see 
the proof of Theorem 17.31 and Remark l7T6l in Section|7] The "interior" Morrey space L''''^ 
is defined to be the set of all functions / e LP{Q.) with finite norm 

IImII^.,, = sup ir-' f lul" 

B,(xo)cn \ JbAxo) 

Moreover, instead of the system ( 11.11 ). one may consider the following system: 

V X (a(jc)V X h) - V(b{x)V ■«)=/ + V x F + V§ in Q.. 
One can show that weak solutions u of the above system are Holder continuous in O if 

g ^6/5.6(i+2W5^ ^g^2,(i+2W2^ ^nd g E L^'™'' for some 5 6 (0, 1). 
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In particular, if / € U'I~(Q), F € U'{Q), and g e L*(Q) for q > 3, then weak solutions 
u € of the above system are Holder continuous in O. Moreover, in that case, we 

have the estimate 

r-'luUB,, + \u\o;B,, < A^('-"'^'ll«ll2;B,. + r^-"'\\f%/2:B,- + r'-"'\\F%,B,. + r'-^'^glU-.B,) , 
whenever B,. - Br{xo) cc Q., where a = a(v, q) e (0, 1) and = A^(v, q). 

Remark 3.8. In Theorem l3.4l one may wish to consider the following problem with non- 
zero Dirichlet boundary data, instead of the problem ( 11.51 ): 

r V X (a(x)V X h) = / + V X g in 
(3.9) I V ■ H = /! in Q, 

\ u - t// on do., 

where one needs to assume the compatibility condition J^h - j^^^ i/f ■ « instead of the 
condition J^h = in Theorem 13.41 If i/f is the trace of a Sobolev function w e W' ''(Q) 
with q > 3, then v := u - w is a solution of the problem (11.51 ) with g and h replaced 
respectively by g and Ii, where 

g := g - aV xw, Ji:^ h-V -w e U'(Q.). 

Notice that Jo^h - 0. Therefore, by the estimate ( 13.51 1 and Morrey's inequality, we have the 
following estimate the weak solution u of the problem i3.9i : 

Mc-(n) ^ N [\\f\\L',n-(n) + Mma) + IMum) + Ikl^i-nn)) , 

where a - a(v, q, Q) e (0, 1) and - N(v, q, Q) > 0. Recall that Q c R-' is assumed to be 
a bounded Lipschitz domain. It is known that if belongs to the Besov space j^^^(c?Q), 
then it can be extended to a function w in the Sobolev space W''*(f2) in such a way that the 
following estimate holds: 

llwlk'^nn) < N\\>/f\\Bl^i^(dn), 
where = N{Q., q); see e.g., Jerison and Kenig lfT3l Theorem 3.1]. Therefore, the follow- 
ing estimate is available for the weak solution u of the problem ( 13. 9t : 

< A^(ll/lb/=(n) + + l|/.||L,(n) + ll<Alb^,„(.a,) , 

where a = a{v,q,D.) e (0, 1) and - N{v,q,Q.) > 0. The above estimate provides, in 
particular, the global bounds for the weak solution u of the problem ( 13. 9t in Q. It seems to 
us that Theorem l3.4l is the first result establishing the global boundedness of weak solutions 
of the Dirichlet problem ( 13. 9t in Lipschitz domains. 

4. Proofs of Main theorems 

4.1. Proof of Theorem 13.11 We shall make the qualitative assumption that the weak so- 
lution u is smooth in O. This can be achieved by assuming coefficient a{x) and the inho- 
mogeneous term / are smooth in O and adopting the standard approximation argument. It 
should be clear from the proof that the constant a and A^ will not depend on these extra 
smoothness assumption. By a standard computation (see e.g., lfT4l Lemma 4.4]), we can 
derive the following Caccioppoli's inequality for u: 

Lemma 4.1 (Caccioppoli's inequality). With u, f, and R as in the theorem, we have 
r IV X H|2 + IV ■ H|2 < A^(r-2 r |„|2 + 11/11^,,,^^^ J ; r = RI3. 
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We take the divergence in the system (11.11 ) to get 

-Ai/r = in Q; i/r • u. 

Denote B(x) = 1 /b{x) and observe that u satisfies 

(4.2) V ■ H = Bi/r in Q. 

Next, we split u=v + w in Br = Br{xo), where r = R/3 and v is a solution of the problem 

V ■ V = BtA - (Bifr),„r in B„ 
V = on dB,, 

where we used the notation 

Jb,(.vo) 

We assume that the function v is chosen so that following estimate, which is originally due 
to Bogovskii UJ, holds for v (see Galdi HI §111.3]): 

(4.3) WmmBr) < N\m - WWIb(B.) < NimWiriB.), Vp e (l, cx,); = Nip). 
Since is a harmonic function, the mean value property of iff yields 

(4.4) m\mB,) < Nr''"-"^\mLHB,,), Vp€(0,H; N ^ N(p). 

Combining the estimates (14.31 1 and (14.41) . and then using ( 14.21) followed by Lemma l^Tl and 
Holder's inequality, we get 

(4.5) I|Vv||l,(b,-) < Nr''''-"^\\u\y(B„.) + Nr\\fh^B,.y, N = A^(v, q). 

By Sobolev inequality, ( 14.31 ), ( 14.21 ). Lemma 1431 and Holder's inequality, we also estimate 

(4.6) ||v||z..(B,.) < Nr\\Vv\y,B,.) < N (HhIL^cb,,.) + r''--"''\\f\y^B,,^) . 

On the other hand, note that w = « - v is a weak solution of the problem 

V X (a(x)V X H") = Vi/r - V X ia{x)W x v) in B,-, 
V • H- = (BiA)x„,,- in B,.. 

We remark that in the proof of | fT4l Theorem 2. 1], we used the condition V ■ m = only to 
establish the following equality (recall the identity ( 11.21 ) above), 

V X (V X h) = -Ah, 

which can be also obtained by merely assuming that V ■ m is constant. Therefore, by lfT4l 
Theorem 2.1 and Remark 2.10], we have (via a standard scaling argument) 

(4.7) r"[wh,B,, < N {r-"^\\wh2^B,.) + r-'^'^ni^'^^B,.-, + r'-^/^'||Vv||w(B,,) , 

where a = a{v, q) e (0,1 - 3/q] and - N{v, q). 

We estimate the RHS of ( 14.71 ) as follows. By the estimate ( 14.6b . we have 

(4.8) \\w\y(B,:) < WuWlHba + \\v\y(B,:) < A^(ll«lb(B3,) + '-'^"'''ll/lb(Bv)) ■ 

By a gradient estimate for harmonic functions followed by ( 14.21 ). Lemma l4n and Holder's 
inequality, we get 

(4.9) ||V^||i,/.(B,.) < Nr'l"-'l-\W\\mB2.-) ^ N[r'^l"-'l^\\u\y^B,,) + r^VW^B,^;) . 
By combining ( 14.7b - ( |4.9b . and ( 14.51 ). we obtain 

r'^lw^B,, < N{r-"%h2^B,,^ + r-''''\\f\\L^iB,.)) ■ 
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By Morrey's inequality followed by ( 14.5b . we also get 

[vVb.. < N[r-'l^-^\\u\y^B^^) + r||/||i,(B,„,) ; = 1 - 3/q. 
By combining the above two estimates and noting that a < ^ - 1-3/(7, we conclude 

(4. 10) r^'iulr.B,, < N [r^^HuWLHB,,) + r--''W\\LHB„.)) ■ 

From the above estimate ( 14.101 ). we can estimate |h|o;B,/4 as follows. For all y € Br/4, the 
triangle inequality yields 

\u{y)\ < \u{x)\ + [u],,Bjr/2r, Vx e 8^/4- 

Taking the average over Br/4 in x, and then using Holder's inequality and ( 14.101 1. we get 

\u(y)\ < |H|2j ' + A^(r-^'/2||H|L2(B„.) + >^-''"\\f\\L'i(B,A) ■ 
Since the above estimate is uniform in y e Br/4, we thus have 

(4.11) |H|o;B,,, < A^(r-'^'ll«lb(B3,.) + r-"'\\f\\LHB,..)) . 

Recall that r = R/3. Therefore, the desired estimate ( 13.2b follows from ( 14.101 ) and ( 14.1 lb 
and the standard covering argument. The theorem is proved. ■ 

4.2. Proof of Theorein l3.4[ We shall again make the qualitative assumption that the coef- 
ficient a{x), the inhomogeneous terms /, g, h, and the domain Q are smooth. By a standard 
elliptic regularity theory, we may then assume that u is also smooth in Q. In this proof, we 
denote by a constant that depends only on v, q, and Q., unless explicitly otherwise stated. 
It should be emphasized that those constants employed in various estimates below, do 
not inherit any information from the extra smoothness assumption imposed on Q; its de- 
pendence on Q. will be only that on the Lipschitz character M, kq of 5f2 and diam Q.. Let us 
recall the following lemma, the proof of which can be found in ITSl . 

Lemma 4.12. Let f e D(Q.), where Q is a domain in R^. Then, there exists F e C°°(Q) 
such that V X F = f in Q.. Moreover, for any p € (I, oo), we have 

llVf |b.(n, < NWfhnny, N = N{p). 

By using the above lemma, we may write / - V x F, where F e C^iQ.) satisfies the 
following estimate 

(4.13) l|VF|b,/2(n) < A^||/lb/^(n); N ^ N{q). 
Notice that u then satisfies 

(4.14) Vx{a{x)Vxu-F- g) = in Q. 
Let (fihe a solution of the Neumann problem 



(4.15) 



A</) = V ■ (a(x)V xu-F - g) in Q, 
d(p/dn = -(F + g) ■ n on 5Q, 



where n denotes the outward unit normal vector of dO.. Recall that ip is unique up to an 
additive constant. We shall hereafter fix (f by assuming j-^ip - 0. 

Lemma 4.16. With u and given as above, we have 

(4.17) V(p^a(x)Vxu-F-g in 
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Proof. First we claim that the boundary condition m = on dQ. implies that 
(4.18) (V xh) ■ « = on dO.. 

To see this, take any surface >S c dO. with a smooth boundary dS c dO.. By Stokes' 
theorem, we then have 

rr (Vxu)ndS =1 udr^O. 

JJS JdS 

Since S is arbitrary and (V x m) ■ n is continuous, we have (V x m) ■ n = on (9fi as claimed. 
Next, we set 

G^V(f- a{x)V xu+F + g. 

The lemma will follow if we prove that G = in Q.. By ( I4.14l i we have V X G = in Q, 
and thus by the condition ( 13. 3t , there exists a potential if/ such that G = Wiff in Q. Then by 
(|4.15t and ( |4.18t , we find that if/ satisfies Aif/ = in O and difjjdn = on 5Q. Therefore, 
we must have G = Vf/r = in Q. The lemma is proved. ■ 



Hereafter, we shall denote A{x) - \fa{x). It follows from (12.11 ) that 

v<A(x)<v"', Vxen. 

Observe that from ( |4.17t we have 

= V ■ (V X m) = V ■ [A(x)(yip + F + g)], 

and thus by Lemma l4.16l we find that satisfies the following conormal problem: 

/ A\v{A{x)V^) = - div {AF + Ag) in Q, 
^ ' \ (A(x)Vip) ■ n = -{AF + Ag) ■ n on dQ. 

In the variational formulation, ( |4.19l l means that we have the identity 

(4.20) 1 AVif ■ = - 1 {AF+ Ag) ■ , e W''^(Q). 
In particular, by using ^p itself as a test function, we get 

\m\LHn) < N lb(n) + ll^llL^cn)) ; = N(v). 
By Poincare's inequality (recall j'o^f - 0) and Holder's inequality, we then have 
||<^||iyi.2(n) < N{\\F\\mn) + Mmsi)) ■ 

Moreover, one can obtain the following estimate by utilizing ( |4.20t and adjusting, for 
example, the proof of dJ] Theorem 8.29] (see and also l22] §VI.10]): 

(4.21) [<p],.,n < N mima) + Wghna^) ; M = Mv, q, Q) 6 (0, 1). 

Then, by Campanato's integral characterization of Holder continuous functions (see 
e.g., 19] Theorem 1.2, p. 70]), we derive from (|4.21| i that 



-i-o) 



(4.22) I \<f-^,,,f <Nr'^^''{\\F\\ 
Jn,(xo) Jci,i 

From the identity ( 14.20b . we also obtain the following Caccioppoli's inequality: 

(4.23) r \V<f\^<Nr-^f |^-^.v„.,f +A^r^-'^*(ll^'lli,(n, + ll^llL,(n))■ 
Ja/2(A■o) -'n,(xo) 
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Setting y = min(ju, l-3/q), and combining (I4.22l i and (14.23b . we get the following Morrey- 
Campanato type estimate for V^: 

(4.24) r |V^|-<A^r'+2y(||^||^^^^^^^||^||^^^^^^)2^ Vxo € Vr € (0, diamQ). 

Having the estimate ( 14.24b together with the boundary condition m = on 5Q, which is 
assumed to be locally Lipschitz, we now derive a global Holder estimate for u as follows. 
Since V ■ « = /i, by ( 11.2b and (14.17b we see that u satisfies 

-Ah = V X (AV^) + V X (AF + Ag) - Vh in Q.. 

By Holder's inequality, we find that (recall y < 1 - 3/^) 

r \F + g\- KNr'^"-^ m\msi) + \\s\\mn)f , V.xo e Q, Vr e (0,diamQ), 
Jn,(A|,) 

where we used the assumption that diamfi < oo. Similarly, Holder's inequality yields 
r \h\^ <Nr^^^~^\\h\\l.,^^^, Vxo€Q, Vr e (0,diamQ). 

Jn,(.>:o) 

Setting G :- A(V(/j + F + g), we find that u satisfies 

(4.25) |-A„=VXG-V. mQ, 

1^ H = on ail, 

where G and h satisfies the following estimate for all jfo e Q and < r < diam Q: 

(4.26) r |G|2 + < Nr'^^y {WWrnn) + WgWmn) + \\h\\mn)f ■ 
Jn,(.vo) 

Observe that the identity ( 14. 17b implies that V x « enjoys the Morrey-Campanato type 
estimate ( 14.24b . The following lemma asserts that in fact, the "full gradient" Vh satisfies a 
similar estimate. 

Lemma 4.27. With u given as above, there exists a — a{v, q, Q) € (0, 1) such that for all 
xq e D. and < r < diam Q, we have 

(4.28) r \Vuf <Nr'^-"{\\Fh,^a) + Mmn) + m\L^n)f■ 
JnAxo) 

Proof. We decompose u = v + w in Qr(xo), where v is the solution of 

-Av = in Q.r{xo), 
V = H on (9Q,(xo). 

Notice that each v' (i = 1,2,3) is a harmonic function vanishing on (dQ.)r{xo) c 5Q. 
By a well-known boundary Holder regularity theory for harmonic functions in Lipschitz 
domains, there exists /3 - y6(Q) e (0, 1) and - N(D.) such that 

(4.29) r \Vv\^ <n(-)^^^^ f |Vv|-, Vp€(0,r]. 
Jnp(xo) ^ r I Jn,(Ao) 

On the other hand, observe that w = « - v is a weak solution of the problem 

-Aw = V X G - V/j in n,.(xo), 
H" = on 5Qr(xo). 

By using w itself as a test function in the above equations and utilizing ( 14.261 ). we derive 

(4.30) r \Vw\^ <N { |G|2 + < A^r'+^r + \\g\\^,^^^^ + Whh^^n)? ■ 
Jn,(A-„) Jn,(xo) 
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By combining ( 14.29b and ( 14.301 1, we get for any p < r, 

r |Vh|2 < i^T^^ r |Vh|- + Nr'^-^ lb(n) + WgWmn) + IMmn)? ■ 

Take any a > such that a < min(j6, y) and applying a well-known iteration argument (see 
e.g. ||9j Lemma 2.1, p. 86]), for all xq e Q. and < r < R < diamfi, we have 

The lemma follows from the above estimate (take R - diam Q.) and the estimate 

(4.31) r |Vh|2 <N f \Gf + |V/2|2 < A^(||F|b(n) + ||g||z.,(n) + Plb(n))' , 
Jn Jn 

which is obtained by using u itself as a test function in ( 14.251 ) and then applying ( 14.26b with 
r = diam Q.. The lemma is proved. ■ 

We now estimate [u]a-Q as follows. Denote by u the extension of u by zero on \ O. 
Notice that u e W^'^(B?) and Vii = ;t'nVH. Then by Poincare's inequality and ( 14.281 ), we 
find that for all x 6 Q and < r < diam Q, we have 

r \u- H,,f < Nr'^-" (IIFIlLnn) + Mmm + m\LHn)f . 

By a standard argument in the boundary regularity theory, it is readily seen that the above 
estimate is valid for all x e Br{xo) and r < 2R, where xq e Q and R - diamQ. Therefore, 
by the Campanato's integral characterization of Holder continuous functions, we find that 
M is uniformly Holder continuous in Br(xo) d f2 with the estimate 

(4.32) [a]aMxo) < N IlLnn, + Mmn) + IMmn)) ■ 
The above estimate ( 14.32b clearly implies that 

(4.33) [H]„;n < m\mn) + Mmn) + IMmn)) ■ 

Finally, we estimate of |M|o;n similar to ( 14.111 ). For xo € Q, the triangle inequality yields 

|h(xo)| < |h(x)| + [Hlccafg^^^.^))/?", VxeQ; = diamQ. 
Taking the average over Q in x, and then using Holder's inequality and ( 14.321 ), we have 

|«(xo)| < (£|Hpj + A^(diamn)" Hi^n) + llglb(n) + IMmn)) ■ 
On the other hand, by (14.31b and the Poincare's inequality, we have 

\u\^<N m~ <Ni\\F\\mn) + \\g\\mn) + \\h\\ 
Jn Jn 

Therefore, by combining the above two inequalities, we obtain 

(4.34) iHlo^n < m\mn> + Mmn) + l|/i|b(n)) ■ 

The desired estimate ( 13.51 ) now follows from ( 14. 33b , ( 14. 34b , ( 14. 13b , and the Sobolev's in- 
equality. The proof is complete. ■ 
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5. Applications 

5.1. Quasilinear system. As a first application, we consider the quasilinear system, 
(5.1) V X m)V X m) - V(S(x, m)V ■ m) = / in O. 

Here we assume y[, S : Q x R^* — > R satisfy the following conditions: 

i) V < S < V"' for some constants v € (0, 1]. 

ii) m and S are Holder continuous in O x R^; i.e. e ^(Q x R^) for ju € (0, 1). 

Theorem 5.2. Let ^ and S satisfy the above conditions and let u 6 F''^(Q) be a weak 
solution of the system (15. Il l with f e U'iQ) for q > 3. Then, we have u G C'-^CQ), where 
a - min(yu, 1 -3/q). 

Proof. By Theorem 13. II we know u e C'^{Q.) for some /3 € (0, 1). Then the coefficients 
a(x) := ^{x,u{x)) and b(x) :- !B{x,u(x)) are Holder continuous with some exponent 
7 G (0, 1). The rest of proof relies on the well-known "freezing coefficients" method in 
Schauder theory and is omitted; c.f. lfT4l Theorem 2.2]. ■ 

Remark 5.3. In Theorem |5^ if one assumes instead that J{,!B e C'^'^'iQ. x M?) and / G 
C*"''''(f2) with keZ+ and fie (0, 1), then one can show that u G C'^^''''(Q); in particular, u 
becomes a classical solution of the system (15.1b . 

5.2. Maxwell's system in quasi-static electromagnetic fields with temperature effect. 

As mentioned in the introduction, the problem (11.5b arises from the Maxwell's system in a 
quasi-static electromagnetic field. Especially, if the electric conductivity strongly depends 
on the temperature, then by taking the temperature effect into consideration the classical 
Maxwell system in a quasi-static electromagnetic field reduces to the following mathemat- 
ical model (see Yin E4i ): 

( H, + Vx (p(m)V xH)^Q, 

[ u, - Au = p(m) |V X H\^, 

where H and u represents, respectively, the strength of the magnetic field and temperature 
while p(m) denotes the electrical resistivity of the material, which is assumed to be bounded 
below and above by some positive constants; i.e., 

(5.4) v<p<v"' for some vg(0, 1]. 

We are thus lead to consider the following Dirichlet problem in the steady-state case: 

V X (p(m)V xH)^0 in Q, 
V ■ H = in Q, 

(5.5) \ H = T on 5Q, 

-Am = p(m) |V X H\- in Q, 

M = on do., 

where we assume that *P and (p are functions in W''''(Q) for ^ > 3. Existence of a pair 
of weak solutions (H, u) was proved in Yin Ii24l and local Holder continuity of the pair 
(//, u) in Q was proved by the authors in lfT4l . Here, we prove that the pair (H, u) is indeed 
uniformly Holder continuous in O. 



Theorem 5.6. Let Q satisfy the hypothesis of Theorem \3.4\ and p satisfy the condition 
(15.4b . Let (H, u) be the weak solution of the problem ( 15.51 ). Then we have (H, u) G C"(Q) 
for some a G (0, 1). In particular, H and u are bounded in Q. 
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Proof. By Theorem [331 and Remark [SJl we find that H e C"{D) for some a e (0, 1) and 
satisfies the estimate 

(5.7) ll^^llc"(n) ^ N\m\w^-nay 

Also, notice from ( 14.28b and Remark [378] that for all xq € Q and < r < diam Q., we have 

(5.8) r |VH|2 < Nr'^^m%, 

Jn,(.vo) 

On the other hand, using the vector calculus identity, 

V ■ (f X G) = (V X f ) ■ G - F ■ (V X G), 
together with the first equation V x (p(m)V x //) = in ( 15.51 ). we find that u satisfies 

-Am = V ■ X (p(m)V X H)) in Q. 
By ( |5.7t and ( 15. 8t , we see that O := // x (piuW x //) satisfies the following estimate: 

(5.9) r |0|- < A^r'+2.||^||4^_ Vxoea Vr € (0, diamQ). 

Jn,-(.v„) 

Therefore, m is a solution of the Dirichlet problem 

-An = V ■ O in Q, 

M = on c)Q, 

where 3> satisfies the Morrey-Campanato type estimate ( 15.91 ) and € W^'''(Q.), and thus by 
a well-known elliptic regularity theory, we have 

In particular, we see that H and u are bounded in Q.. The proof is complete. ■ 

Remark 5.10. In Theorem |531 if one assumes further that p 6 C'^(R), where k G Z+, then 
by Theorem 15.21 and the bootstrapping method, one finds that H e C^'"'(Q) n C"(Q) and 
u e C*+' "(Q) n C"(Q); see IH Theorem 3.2 and Remark 3.3]. In particular, if p 6 C^(R.), 
then the pair {H, u) becomes a classical solution of the problem (15.51 1. 

6. Green' s FUNCTION 

In this section, we will discuss the Green's functions (more appropriately, it should be 
called Green's matrices) of the operator L in arbitrary domains. Let Z be any subset of Q 
and u be a y'''(Q) function. Then we shall say u vanishes on S (in the sense of Y^'^iQ.)) if 
M is a limit in Y^-^{Q) of a sequence of functions in C^(f2 \ S). 

Definition 6.1. We say that a 3 x 3 matrix valued function G{x,y), with entries Gij{x,y) 
defined on the set [{x,y) e Q. x Q. : x y], \s ?l Green's function of L in f2 if it satisfies the 
following properties: 

i) G{-,y) € W];l{a) and LG{-,y) = 6yl for all y e Q, in the sense that for = 1, 2, 3, 

r fl(V X G{;y)ek) ■ (V x ^) + b{V ■ G(-,y)ct)(V ■ <l>) = /(y), V0 € C^(Q), 
Jn 

where Cjt denotes the k-\h unit column vector; i.e., ei - (1,0, 0)^, etc. 

ii) G{-,y) € F''2(Q \ Br(y)) for ally e Q and r > and G(-,y) vanishes on dQ. 
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iii) For any / € CJ^(il), the function u given by 

u{x):= f G(y,x)f(y)dy 
Jn 

is a weak solution Yq'^(D.) of the problem i2AV . i.e., u belongs to Yq^(Q.) and satisfies 
Lu - f in the sense of the identity ( 12.3b . 

We note that part iii) of the above definition gives the uniqueness of a Green's matrix; 
see Hofmann and Kim jlSI- We shall hereafter say that G(x,y) is the Green's matrix of L 
in O if it satisfies all the above properties. Then, by using Theorem [TT] and following the 
proof of lfT2l Theorem 4.1], we obtain the following theorem, where we use the notation 

a A b :- min(a, b), aV b :- max(fl, b), where a,b eR. 

Theorem 6.2. Let Q. be a (possibly unbounded) domain in R^. Denote d^ :— dist(x, 5Q) 
for X € Q; we set dx — co ifO. — M?. Then, there exists a unique Green's function G{x,y) 
of the operator L in Q, and for all x, y € Q satisfying Q <\x — y\ < d^ A dy, we have 

(6.3) \G(x,y)\ < N\x - y\-\ where N ^ N{v) > 0. 

Also, we have G(x, y) — G(y, x)^ for all x,y e Q. with x + y. Moreover, G(-, y) 6 C"'(Q \ {y\) 
for some a — a(v) e (0, 1) and satisfies the following estimate: 

(6.4) \G{x,y) - G{x',y)\ < N\x - x'\"\x - y\-^-", where N = A^(v) > 0, 
provided that \x — x'\ < \x — y\/2 and \x —y\< dx A dy. 

Next, we consider the Green's functions of the system (11.3b . 

Definition 6.5. We say that a 3 x 3 matrix valued function G{x, y), which is defined on the 
set [(x,y) e Q X Q : X 9^ y), is a Green's function of the system (11.3b in Q if it satisfies the 
following properties; 

i) G(- , y) € (Q) for all y € O and for = 1 , 2, 3, we have 

r fl(VxG(-,y)ei)-(Vx^) = /(3;), V0 e C^(Q), 
Jn 

r G(-,3;)e,-V<A = 0, V«A e Co"(Q), 
Jn 

where Ck denotes the k-ih unit column vector; i.e., e\ - (1,0, 0)^, etc. 

ii) G(-,y) e Y^--(Q. \ B,{y)) for all y 6 Q and r > and G{-,y) vanishes on dQ.. 

iii) For any / e 0{Q), the function u given by 

H(x):= r G(y,x)f(y)dy 
Jn 

is a weak solution in Y^'^{D.) of the problem 

( V X (fl(x)V xu)^f in n, 
i V ■ H = in 

y u - on dD., 

1 2 

that is, u belongs to Y^' (Q) and satisfies the above system in the sense of the identities 
(l23T l and with g = and /j = 0. 
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Then by the same reasoning as above, Theorem l6.2l also applies to the Green's functions 
of the system ( 11.31 ). Moreover, in the case when Q is a bounded Lipschitz domain satisfying 
the condition ( 13.31 ), a global version of estimate ( 16.3b is available thanks to Theorem 13. 41 
and lfT6l Theorem 3.13]. 

Theorem 6.6. The statement ofTheorem \6.2\ remams valid for the Green 's functions of the 
system ( 11.3b . Moreover, if we assume that Q is a bounded Lipschitz domain satisfying the 
condition (13.3b . then for all x, y 6 f2 with x + y,we have 

\G{x,y)\ < N[d, A \x-y\]''[dy A \x - y\f\x - y\-'-^\ 

where a - a(v, Q.) e (0, 1) and N - N(v, Q). 

Remark 6.7. Theorem l6.2l in particular establishes the existence of the Green's function 
of the operator L in R^, which is usually referred to as the fundamental solution of the 
operator L. Notice that in that case, we have the pointwise estimate ( 16.3b available for all 
x,y e R-^ with x y, and estimate ( 16.41 ) for all x, x' satisfying \x-x'\ < \x-y\/2. The various 
estimates for the Green's function that appears in lfT2l Theorem 4.1] are also available in 
Theorem l6.2l 

7. Associated parabolic system 
In this separate and independent section, we consider the system of equations 

(7.1) H, + V X (fl(x)V X h) - V(/7(x)V ■ h) = / in Q X (0, T), 

and prove that weak solutions of the system ( 17.11 ) are Holder continuous in Q x (0, T) pro- 
vided that / satisfies some suitable condition, which is an extension of ifTTl Theorem 3.1], 
where it is shown that weak solutions of the following system are Holder continuous; 

(7.2) H, + V X ia{x)V X h) = 0, V ■ « = in O x (0, T). 

As mentioned in the introduction, the above system arises naturally from Maxwell's equa- 
tions in a quasi-static electromagnetic field. More precisely, let cr(x) denote the electrical 
conductivity of a material and the vector H{x, t) represent the magnetic field. It is shown 
in Landau et al. 1(211 Ch. VII] that in the quasi-static electromagnetic fields, H satisfies the 
equations 

+ V X (;i:V X /y) = 0, V ■ H = in O X (0, T), 

which is a special case of the system ( 17.1b . Also, in this section we study the Green's 
functions of the system ( 17.1b and the system ( 17.21 ). by using recent results from ||2l|3l. 

7.1. Notation and definitions. In this section, we abandon some notations introduced in 
Section |3] Instead, we follow the notations of Ladyzhenskaya et al. |[T9l with a slight 
variation. We denote by Qj the cylindrical domain Q x (0, T), where T > is a fixed 
but arbitrary number, and St the lateral surface of Qt; i.e., St - <5Q x [0, T}. Parabolic 
function spaces such as L^AQt), LgiQj), W'/{Qt), W'/{Qt), V2(Qt), and V^-\Qt) are 
exactly those defined in Ladyzhenskaya et al. |fT9l . We define the parabolic distance 
between the points X - {x, t) and Y - {y, s) by 

\X-Y\p :=max(|x-3;|, ^7^) 

and define the parabolic Holder norm as follows: 

\u{X) - u{Y)\ 

\u\al2.a:Q = [u]a/2,a-Q + \u\q;Q := SUp — — — — + SUp|M(X)|. 

XJeQ |A - -Tip XeQ 
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We write Vm for the spatial gradient of u and m, for its time derivative. We define 

Q;(X) = X (r - r\ t), Q,iX) = B^ix) x{t-r\t + r\ 

We denote by ^ the operator d, + L; i.e., 

ifu := H, + Lh = H, + V X (fl(x)V x «) - V(/7(x)V ■ «), 

and by '.if the adjoint operator -d, + L. For a cylinder Q of the form Q x (a, h), where 

1 



-oo < fl < ^ < oo, we say that m is a weak solution in y2(0 (Vj'^Q)) of = / if 
« e V2(0 (Vi'^CQ)) and satisfies the identity 



- r H-0,+ r fl(vxH) (vx^)+Mv-H)(v-^)= r z-^, m^^c^^q). 

Jq Jq Jq 

Similarly, we say that « is a weak solution in VoSQ) (Vj'^CG)) of '^m = / if m e V^SQ) 
(V,'"(g)) and satisfies the identity 

r u ■(/>,+ [ fl(VxH) (Vx0) + /7(V H)(V -^)= [f (f,, V(f,eC^iQ). 

Jq Jq Jq 

7.2. Holder continuity estimates. The following theorem is a parabolic analogue of The- 
orem [TT] However, it should be clearly understood that in the theorem below, the coeffi- 
cients a and b of the system ( 17.11 ) are assumed to be time-independent. 



Theorem 7.3. Let Qt = Q x (0, T), where Q. be a domain in . Assume that a(x) and 
b(x) are measurable functions on Q satisfying ( 12. 11 ). Let u be a weak solution in V2(Qt) of 
the system ( 17.11 ) with f e LqiQj) for some q > 5/2. Then u is Holder continuous in Qt, 
and for any — Q'j^iXo) CC Qt, we have the following estimate for u in 2^^2-' 

(7.4) R"[uU„/2,Qj^^^_ + |H|o;e,,, < A^(^"'^'ll«lk(Qp +^'"''*II/IIl,(Gp)> 

where a - a{v, q) e (0, 1) and N - N(v, q) > 0. 

The proof of the above theorem will be given in g7.4l below. 

Remark 7.5. As in IfTTl Theorem 3.2], one can consider the case when the coefficients a 
and b of the system ( 17. Il l are time-dependent but still have some regularity in f-variable. 
For a measurable function / — f(X) = f{x, f), we set 

w^C/) := sup sup ] [ [ \f(y,s)- f^riy)\dyds, V5 > 0, 

X=(t,x)ER* r<6 \Qi\^)\ Jt-r- Jb,(x} 

where fj{y) - j^^, f(y,s)ds. We say that / belongs to VMO, if limg^o ajs(f) - 0. 
Assume that the coefficients a{x, t) and b{x, t) are defined in the entire space and belong 
to VMOr. Let H e V2{Qt) be a weak solution of the system 

H, + V X (a{x, r)V X h) - V(/7(x, f)V ■ u) = f in Qt, 

where / € Lq{QT) with q > 5/2. Then one can show that u is Holder continuous in Qt. The 
proof is very similar to that of ifTTl Theorem 3.2]. Also, as is mentioned in Remark l376l one 
may assume that a is a 3 x3 (possibly non-symmetric) matrix valued function satisfying the 
uniform ellipticity and boundedness condition; see |i7j and also consult IfTSl for treatment 
of non-symmetric coefficients. 



ELLIPTIC SYSTEMS WITH MEASURABLE COEFFICIENTS 



17 



Remark 7.6. In Theorem 17. 3 1 instead of assuming that / € LqiQj), one may assume that 
/ belongs to the mixed norm space Lq^riQi) with suitable q and r. In fact, one may assume 
that / belongs to the Morrey space, m^'^I'^'^^^^^^^I'^ with 6 € (0, 1), where M'' '? is the set of 
all functions / € Lp{Qj) with finite norm (c.f. Lieberman ||22] §VI.7]) 

\ Jqaxo) I 



\\U\\mp.i = sup 



Then, instead of the estimate ( 17.22b in the proof of Theorem l7.3l we would have 

IVh-i^ < A^II/llLMer(x)) ^ a^'-'"''II/II'm.»/-o,3.»/,- 



I 



>QAX) 

The rest of proof remains essentially the same. 

7.3. Green's function. Let f/ = Q x R be an infinite cylinder in with the base Q being a 
(possibly unbounded) domain in and let dU be its (parabolic) boundary dQ. x R. Let 
S cz Q and m be a W2^(Q) function. We say that u vanishes (or write u = 0) on >S if m is a 
limit in W2^(Q) of a sequence of functions in C^iQ \ S). 

Definition 7.7. We say that a 3 x 3 matrix valued function G(X, Y) - G(x, t, y, s), with 
entries Gij{X, Y) defined on the set [(X, Y) e U x U : X Y], is a Green's function of the 
operator ^ in [/ if it satisfies the following properties: 

i) G(-, Y) 6 Wl'l^iU) and ^G(-, Y) = dyl for all Y eU,m the sense that for /t = 1, 2, 3, 
the following identity holds for all (f>eC^{U): 

[ -G{; Y)ek ■ 4>, + fl(V X G{;y)ek) -(7x0) + ■ G(-,y)e,)(V ■ 0) = /(F), 
Ju 

where Ck denotes the A;-th unit column vector; i.e., ei - (1,0, 0)^, etc. 

ii) G(-, Y) € Vl'^iU \ Qr{Y)) for all F € and r > and G(-, Y) vanishes on dU. 

iii) For any / e C^{U), the function u given by 

H(X):= r G(Y,X)f(Y)dY 
Ju 

is a weak solution in y2'°(t/) of 'Ji'u - f and vanishes on dU. 

We note that part iii) of the above definition gives the uniqueness of a Green's function; 
see 121. We shall thus say that G{X, Y) is the Green's function of ^ in t/ if it satisfies the 
above properties. By Theorem |7.3| and ^ Theorem 2.7], we have the following theorem: 

Tlieorem 7.8. Let U — Q. x W be an infinite cylinder, where the base Q. is a (possibly 
unbounded) domain in M?. Then the Green 's function G{X, Y) ofj^ exists in U and satisfies 

(7.9) G{x,t,y,s)^G(x,t-s,y,Oy, Gix,t,y,Q) = for t < 0. 
For all f e C^(C/), the function u given by 

(7.10) H(X):= r G(X,Y)f{Y)dY 

Ju 

is a weak solution in vY'iU) of — f and vanishes on dU. Moreover, for all g € L'^(Q), 
the function u(x, t) defined by 

u(x,t) := I G{x,t,y,0)g(y)dy 
Jn 
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is a unique weak solution in Vj^'CGr) of the problen^ 

and if g is continuous at xq & Q. in addition, then we have 

lim u{x, t) - g{xQ). 

(x,r)^(xo,0) 
xeSl, t>0 

Remark 7.11. The identity G(x, t,y, s) - G{x, t - s,y, 0) in Theorem 17 . 81 comes from the 
fact that ^ has time-independent coefficients; see Q. The function Kr(x,y) defined by 

(7.12) K,{x,y)^Gix,t,y,0), x,y e Q, t>0 

is usually called the (Dirichlet) heat kernel of the elliptic operator L in O. It is known that 
K, satisfies the semi-group property 

K,+Ax,y)= f K,{x,z)K,(z,y)dz, \/x,yen, Vr,i>0, 
Jn 

and in particular, if Q = R^, then we also have the following identity: 

K,(x,y)dy = I, Vx6R\ Vf > 0, 

where / denotes the 3x3 identity matrix; see ID Theorem 2. 11 and Remark 2. 12]. 

The following theorem is another consequence of Theorem l7.3t see ^ Theorem 2.11]. 

Theorem 7.13. Let Kt{x,y) be the heat kernel for the operator L in 9? as constructed in 
Theorem \7.8\ Then we have the following Gaussian bound for the heat kernel: 

\K,(x,y)\ < Nt'^'^exp{-K\x~y\^/t}, Vf > 0, x,y e R\ 

where N - N(v) > and k - k{v) > 0. 

Next, we consider the Green's functions of the system ( 17.21 ). 

Definition 7.14. We say that a 3 x 3 matrix valued function G(X, Y) - G(x, t, y, s), with 
entries Gij(X, Y) defined on the set {(X, Y) e U x U : X ^ Y], is a Green's function of the 
system (17.21) in U if it satisfies the following properties: 

i) G(-, Y) e W}f(U) for all F e f/ and for /t = 1, 2, 3, we have 



[ -G{ , Y)et ■ (p, + fl(V X G(;y)e,) -(7x0) = 4>'{Y), V<^ e C^(U), 
Ju 

G{ , Y)et ■ ViA = 0, V(A e C^{U), 

u 

where Ck denotes the A:-th unit column vector; i.e., Ci - (1,0, 0)^, etc. 

ii) G(-, Y) e y^'"(t/ \ QAY)) for all 7 e and r > and G{ , Y) vanishes on dU. 

iii) For any / e C^(U) satisfying V ■ / = in [/, the function u defined by 

H(X):= r G(Y,X)f{Y)dY 
Jn 

is a weak solution in V,' '"([/) of the problem 



X 



-M, 4- V X (fl(x)V xu) = f, V • H = 0, 



'See, Ladyzhenskaya et al. [19] §111.1] 
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that is, u belongs to (U), vanishes on dU, and satisfies the above system in the 
sense of the following identities: 



Jn 

It can be easily seen that existence of the Green's function of the system ( I7.2l i in U 
follows from ITtI Theorem 3.1] and 12] Theorem 2.7], and that it satisfies the relations 
( |7.9l l in Theorem l7.8l We shall say that K, defined by the formula ( 17.12b is the (Dirichlet) 
heat kernel of the elliptic system ( 11.31 ) in Q. Then it satisfies the statement in Remark l7.11l 
as well as that in Theorem 17. 131 If we assume further that f2 is a domain satisfying the 
hypothesis of Theorem 3.5, then we have the following result, which is an easy conse- 
quence of in Theorem 3.6] combined with Theorem l3.4l and Lemma 4.4] (see also [S] 
Remark 3. 10]): 

Theorem 7.15. Let C/ = Q x R with Q. satisfying the hypothesis ofTheorem \3.4\ Then the 
heat kernel K,(x,y) of the system ( 11.31 ) exists in Q. Moreover, for all T > there exists a 
constant N — N{v, f2, T) such that for all x,y e D. and < t < T, we have 



where k — k(v, Q.) > and a — a(v, O) € (0, 1) are constants independent of T, and we 
used the notation a A b — min(fl, b), aW b — max(fl, b), and d^ — dist(x, dQ). 

7 A. Proof of Theorem |731 We follow the strategy used in 1 17 j. As before, we shall make 
the qualitative assumption that the weak solution u is smooth in Qt- Let us first assume 
that / = and consider the homogeneous system 

(7.16) H, + Lh := H, + V X (a(jc)V X m) - V(/7(x)V ■ h) = in gr- 

The proof of the following lemma is very similar to that of IfTTl Lemma 3.1 - 3.3], where 
we strongly used the assumption that coeflicients of the operator are time-independent. 

Lemma 7.17. Let v e V2(2^,.), where Q^^^ — Q^^X^q) with A > I, be a weak solution of 
Vt + Lv — in 2^7r- Then we have the following estimates: 



where N - N{v, A) > 0. 

The proofs of following lemmas are also standard in parabolic theory and shall be omit- 
ted; see e.g., ||2l Lemma 2.4 and 3.1] and also |3] Lemma 8.6]. 

Lemma 7.18. Let u e V2{Q.^), where - Q,^{Xq), be a weak solution ofut + Lu - f in 
Q^. Then we have the estimate 
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where N - N(v) > 0. 

Lemma 7.19. Let u e V'2(G^,.), where Qj^. = 2^,.(^o) with A > \, be a weak solution of 
u, + Lu - f in Q^^. Then we have 



[ \u(;t)\^+ [ m^<N{r-'- [ \u 
Jb,. Jot \ Joi 



sup 

io-r-<i<io -Jb,- JQ] 



ll-^lli.lO/7(Gj,) 



where N - N{v, A) > 0. 



With the above lemmas and Theorem 13.11 at hand, we now proceed as in the proof of 
ifTTl Theorem 3.1] (see also proof of |fT8] Theorem 3.3]) to conclude that any weak solution 
V € V2(Qt) of the system ( 17.16b is Holder continuous in Qt and satisfies the estimate 

(7.20) [vV,/2;Q,„ < NR-''^-^'\\vh,jQ-^y, e« = Q-^(Xo), 

where fi - //(v) e (0, 1) and - N{v) > 0. There is a well-known procedure to obtain 
Holder estimates for weak solutions of the inhomogeneous system u, + Lu - f from 
the above estimate ( 17.20b for weak solutions of the corresponding homogeneous system 
u, + Lu - 0, which we shall demonstrate below for the completeness. For X € Q^i^{Xq) 
and r e (0,/?/4], we split u - v + w m Q^{X), where w is the unique weak solution 
in ''(Gr (^)) of w, + Lw = / in Q'^(X) with zero boundary condition on the parabolic 
boundary dpQ^{X). Then, v - u - w satisfies v, + Lv = in Q^iX), and thus, for < p < r 
(c.f. m Eq. (3.9)]), we have 



(7.21) r |Vm|2<2 r |Vv|-+2 f 

<N(p/rf^^'' f |Vv|2 + 2 r IVh'P 

<N{plrf^^i' r \Vu\^+N r IVh-P. 

jQTiX) Jqax) 

Choose p e (5/2, q) such that a :- 2 — 5 1 p < ^i. By the energy inequality and a parabolic 
embedding theorem (see HH §11.3]), we get (c.f. ^ Eq. (3.10)]) 

•J Qr (X) 

Combining dTTTTb with (|7^ , we get for all p < r < /?/4, 



r \Vu\- <N{plrf^-'' r \Vu\^+Nr^ 

•JQ-piX) JQTiX) 

Then, by a well known iteration argument (see e.g., (|9] Lemma 2.1, p. 86]), we have 
r |Vh|2 < N(r/Rf^^- f |Vh|2 + A^r^+2„||^||2 

Jerm Je„„(X) ^"^^""^ 

By Lemma |7.18l the above estimate, and Holder's inequality, we get 

LJ" - < Nr'-'" (/^-^-^"IIV«llLe,,(X)) + ll/llt(G, J ■ 
Then, by Campanato's characterization of Holder continuous functions, we have 
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By Lemma l7.19l and Holder's inequality (recall a = 2 - 5//?), we then obtain 

(7.23) ^"[«]„,a/2;G,,, < (/?-'/'||h||l,(q,) + ^'-'^''ll/llL,(e,)) • 
Similar to (14.1 11 1, we then also obtain 

(7.24) |H|o;e,-,, < n{r-'''\\u\\,^jq-^, + R'-''V\K,Q„,) . 

Finally, the desired estimate ( 17.4b follows from ( 17.23b . (17.24b . and the standard covering 
argument. The theorem is proved. ■ 

8. Appendix 

8.1. Existence of a unique weak solution of the problem ( 12.4b . We prove existence of a 

1 ^ 

unique weak solution in }'y"'(Q) of a more general problem 
(8.1) 



V X (a(ji)V X h) - V{b(x)V ■u)^f + VxF + Vg in Q, 
M = on do., 

where / € L^'\^l) and F,ge L-(Q.). We say that a function m is a weak solution in Fq"(Q) 
of the problem ( 18.1b if u that belongs to ^^'^(Q) and satisfies the identity 

1 fl(V X h) ■ (V X v) + b{W ■ h)(V ■v)= 1 f-v + F-Vxv + gV-v, Vve Cq{Q.). 
Jn Jn 

Notice that the inequality (12.2b imphes that the biUnear form 
(8.2) <M,v) = {u,v}h = y f Vm' ■ Vv' 

12 T 

defines an inner product on H :— F,,' (Q.) and that H equipped with the above inner 
product is a Hilbert space. We define the bilinear form associated to the operator L as 

B[u, v] := I fl(V X m) ■ (V X v) + ■ m)(V ■ v). 

Then, in light of the identity ( 11.21 ). we find that 

1 |VxH|2 + |V-Hp = 1 |Vh|', VHe//. 
Jn Jn 

It is routine to check that the bilinear form B satisfies the hypothesis of the Lax-Milgram 
Theorem. On the other hand, by the inequality ( 12. 2b . the linear functional 



F(v):= [ / v + F-Vxv + gV v 

Jn 

is bounded on H. Therefore, by the Lax-Milgram Theorem, there exists a unique element 
u e H such that B[u, v] = F(v) for all v e H, which shows that u is a unique weak solution 
in Y^'^iQ.) of the problem dO- ■ 
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8.2. Existence of a unique weak solution of the problem (11.5b . We shall assume that 
/ e H(,/5{D.) and g € L-{D.). First, we consider the case when li - and construct a weak 
solution in Y^"{Q.) of the problem ( 11.51 ) as follows. Let H be the completion of D{Q.) (see 
Section |2] for its definition) in the norm of y''2(0). Then H c Y^^-(n)^ and as above, it 
equipped with the inner product ( 18.21 ) becomes a Hilbert space. We define the bilinear form 
BonHhy 

B[u,v] := f a(V xu) ■ (V x v). 

Then the bilinear form B satisfies the hypothesis of the Lax-Milgram Theorem. We also 
define the linear functional F on H as 

F(v):= r f-v + g-(Vxv). 
Jn 

One can easily check that F is bounded on H. Therefore, by the Lax-Milgram Theorem, 
there exists a unique element u e H such that B[u, v] = F{v) for all v e H. In particular, u 
satisfies identities ( 12.5b and ( 12.61 ) with h - 0. Therefore, « is a weak solution in yjJ'~(Q) of 
the problem ( 11.5b in the case when h = 0. 

Next, we consider the case when h + 0. In this case, we assume further that f2 is a 
bounded Lipschitz domain so that in particular, we have Fq'^(Q) — Wq'^(Q). For h 6 L^(Q) 
^ /; = 0, let V e Wq' (Q) be a solution of the divergence problem 

J V ■ V = /j in O 
I V = on 

that satisfies the following estimate (see e.g., Galdi 18] §111.3]) 

l|Vv|L.(n) < mhWum. N = N(Q). 

Let w be a solution in Yq'(Q.) of the problem ( 11.5b with g - aV x v in place of g and h - 0, 
which can be constructed as above. Then, it is easy to check that u := v + h" is a solution 
in ^^'^(Q) = W„'^(0) of the original problem ( fT3] l. 

Finally, we prove the uniqueness of weak solutions in Yq^(Q.) of the problem (11.5b under 
the assumption that O is a bounded Lipschitz domain. Notice that in that case we have 
^^'"(Q) = Wq'^(Q). Suppose u and v are two weak solutions in Wg'^(O) of the problem 
( 11.5b . Then the difference w = « - v is a weak solution in W,]'^(0) of the problem ( 11.5b 
with f - g - and h - 0. By the identity ( |2.6b . we find that w e H; see e.g., Galdi lH) 
§111.4]. Then by the identity ( 12.51 ). we conclude that w - 0, which proves the uniqueness 
of weak solutions in Fg'^(Q) of the problem ( 11.51 ). ■ 
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